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Abstract. Bragg scattering with linearly polarized light can be used to separately 
measure the density and the spin structure factor of a two spin component atomic gas 
by looking at the dependance of the scattering intensity on the polarization of the laser 
fields. Both stimulated and spontaneous scattering are discussed. Explicit results for 
different spin configurations are given. 

PACS numbers: 67. 75.40. Gb 52.38.Bv 
1. Introduction 

In recent years, a strong activity has been concentrated on the study of atomic gases 
with spin degrees of freedom [Tj. In this context, optical diagnostic techniques appear 
very promising, as they provide the possibility of selectively addressing the different 
spin states by playing with the light frequency, as in the first observation of vortices in 
Bose-Einstein condensates |2J, or with the polarization degrees of freedom, as in recent 
experimental studies of the magnetization in a spin-1 atomic cloud [HUE]. 

Among the many effects under investigation, the superfluid transition of two spin 
component Fermi gases close to a Feshbach resonance [HJ Ej has received a particular 
interest from both the theoretical and the experimental side since atomic gases are 
clean systems where it is possible to quantitatively test different theories for strongly- 
interacting many-fermion systems jTj. One of the key observables for the microscopic 
characterization of such a system are the density-density correlation functions. In 
particular the opposite spin correlation function Sj_(x, x') = (n + (x) n_(x')) carries 
information on the pairing effect which is responsible for superfluidity jH]. Recent 
experiments using molecular spectroscopy techniques have measured its local value 
S_i_(x = x') 0, but no data are yet available at different points r ^ r'. 

Scattering techniques are widely applied to the study of condensed matter systems, 
since they allow for a precise characterization of the structure and the dynamics down 
to length scales of the order of the wavelength used in the experiment. Remarkably, 
the typical microscopic length scale of ultracold atomic gases is often on the order of 
the visible wavelength, and this has allowed the use of stimulated Bragg scattering 
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techniques, e.g. to measure the structure factors (i.e. the Fourier transform of the 
density-density correlation function) in clouds of spinless bosonic atoms [TU] . 

In the present paper, we propose to take advantage of the polarization degrees 
of freedom of light in order to perform Bragg scattering experiments which are able to 
selectively address the total density and the spin density of the atomic gas, and then give 
information on the same- and opposite-spin density- density correlations. After a brief 
review of the light-matter interaction for spin-1/2 atoms, we shall show that a stimulated 
Bragg scattering experiment analogous to the one of [10] is able to measure the response 
function to perturbation operators which are a combination of the total and the spin 
densities with polarization-dependent weights. The fluctuation-dissipation theorem then 
relates the absorbed energy to the structure factor. A direct measurement of the 
structure factor, not relying on the fluctuation-dissipation theorem, can in principle 
be performed by means of a spontaneous Bragg scattering experiment using a single 
incident beam and a polarization-sensitive detection. 

2. The optical response of spin 1/2 atoms 

Consider a gas of spin F g = 1/2 fermionic atoms with an optical transition at frequency 
u eg to an excited state of spin F e . In the simplest case F e = 1/2, the light-matter 
coupling Hamiltonian in the rotating-wave approximation can be written in the compact 
form [1]: 



in terms of the electric dipole moment d of the optical transition. Here, cr % are the Pauli 
matrices [i = {x,y,z}), and Ei are the electric field operators; ^ g , a and \l/ ejQ are the 
atomic field operators for the Zeeman sublevels a, [3 — ± of respectively the ground g 
and the excited e states. In the F e = 3/2 case, the structure is analogous, the index 
a now running over the four Zeeman sublevels a = 3/2, 1/2, —1/2, —3/2 of the excited 
state, and the Pauli matrices being replaced by the ones of the relevant Clebsch-Gordan 
coefficients. 

Assuming that the light interacting with the atoms is far from resonance with the 
optical transition, we can adiabatically eliminate the excited state and write an effective 
Hamiltonian involving the ground state only. If no external field is present which breaks 
the rotational symmetry, the Zeeman sublevels of both the ground and the excited 
states are degenerate, so that the different tensorial components of the light and matter 
fields can be separated, as originally done in [EQ. For our F g = 1/2 case, the effective 
Hamiltonian has then the form |12j : 



where repeated indices are meant to be summed over. The effective coupling constants 
fn tS depend on the electric dipole moment d of the optical transition and on the detuning 




(1) 




dr (/„ 6a <W + i fs $00 *J>) ^ 
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from the excited state. For the simplest case oi a F g = 1/2 — > F e = 1/2 transition, one 
has f n — fs — \d\ 2 /h{u)L — w ej ), while for a F g = 1/2 — > F e = 3/2 transition one has 
f n = —2 f s = 4 \d\ 2 /3h(uL — Meg)- If several transition are involved, their contributions 
to / n>a have to be algebraically summed. 

By means of elementary vectorial analysis, this Hamiltonian can be rewritten as: 

H ^ = \J dr [fn n(r) J(r) + f s S(r) ■ B(r)] (3) 

in terms of the intensity J(r) and the pseudo-magnetic field B(r) defined as [T2"] : 

J(r) = E f (r) • E(r) (4) 
B{v) = zE f (r) xE(r). (5) 

Note how the coupling of the light field to the electric dipole of the far-off resonant 
optical transition results in the effective potential (j3J) which not only couples to the atom 
density h(r) = *| iQ! (r) ^ g , a {r), but also to the spin density Si(r) = a l af3 ^\ t Jr) ^ g ,p{ v )- 
The physical meaning of this coupling can be easily understood at least in some most 
significant cases. For linearly polarized light, the vector product vanishes, so that the 
optical potential couples to the total density only. On the other hand, for a cr + -polarized 
light, the pseudo-magnetic field is B = —Iz. In the case of a transition to a F e = 1/2 
excited state for which f n — f s , it is easy to see from (jHJ) that the optical potential (jHJ) 
created by cx + -polarized light acts on the — state only. No optical transition can in fact 
be driven on atoms in the + state by cr+ light. 

3. Stimulated Bragg scattering 

Consider a stimulated Bragg scattering configuration in which a pair of laser beams of 
wavevectors k 12 are sent onto the atomic cloud with polarization vectors E 12 . If the laser 
frequencies uj\ and U2 = iO\ — Au are well detuned from atomic resonance, spontaneous 
emission can be neglected, and the electric field operators in the Hamiltonian (J3J) 
can be replaced by the classical C-number amplitudes given by the superposition 
E(r) = Ei(r,t) + E 2 (r,t) of the two beams Ei, 2 (r, t) = E? )2 (r) e ikl - 2r e - ^ 1 - 2 *. The spatial 
dependence of E° 2 (r) accounts for the finite waist of the two beams. 

The spin-dependent optical potential acting on the atoms then reads: 

H eS = i J dr [f n (Ef(r) • E°(r) + E°*(r) • E°(r)) n(r) + 

+ if, (ET x E° + E°* x E°) • S(r) + 

+ f n (E°* (r) • E°(r) e -*(*H-Ao,t) + h c ^ - (r) + 

+ i f s (e°* x E° e -i(i-*+ A "t) + h.c) • S(r) . (6) 

The first two terms describe the optical potential created by each laser beam separately, 
and has a negligible effect on the atoms as long as the waist of the beams is much wider 
than the atomic cloud. The third and fourth terms describe the potential due to the 
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interference pattern of the two beams. This has a wavevector q = k 2 — and moves 
at speed v = —Au/q in the direction of q. 

A very interesting case is when both beams have linear polarizations along the 
directions e 12 making a relative angle 9 such that cos# = ei • e 2 . Neglecting the effect 
of the finite waist of the beams, the optical potential Hamiltonian @ can be rewritten 
as: 



g -i(q-r+Au;t) _|_ 



/„ (e 2 ■ ex) n(r) + i f s (e 2 x e : ) • S(r) 

+ H.c. (7) 
which has the standard form of a linear response Hamiltonian: 

# eff = Voie-^' + H.c, (8) 
the perturbation operator A being: 

A = f n cos6»n(q) +if s sm9 S z (q), (9) 

where 

n(q) = f dre" ,q - r n(r) (10) 



S(q)= J rfre~ iq - r S(r) (11) 

are the Fourier components of the total and spin densities. 

Depending on the relative weight of the coupling to the density and to the spin 
density, different families of elementary excitations [HI E3 EJ are addressed by the 
perturbation A. The Anderson-Bogoliubov phonon ^H] couples in fact to the density, 
while pair-breaking excitations couple to the spin density. Here, one can go from one 
case to the other simply by varying the angle 9 between the linear polarizations ei, 2 . 
If ei )2 are parallel (9 = 0), the field is everywhere linearly polarized, with a sinusoidal 
spatial dependance. This means that the vector product in ((7j) vanishes and the optical 
potential couples to the total density only. On the other hand, if e 12 are orthogonal 
9 = 7r/2 (lin _L lin configuration JH]), the field intensity J(r) is spatially constant, but 
the polarization periodically varies in space as cr + — > (x + y) — > er_ — > (x — y) — > <r + 
giving a sinusoidally varying pseudo-magnetic field B(r) polarized along the direction 
perpendicular to e 12 . In the following we shall take this direction as the quantization 
z axis: in this way, the optical potential results always diagonal in the ± basis of the 
atomic states. If a different quantization axis is chosen, e.g. orthogonal to ei x e 2 , the 
coupling of the pseudo-magnetic field £>(r) to the spin would result in off-diagonal terms 
coupling the two different spin states in a Raman-like way. 

It is a general result of linear response theory that the imaginary part 
Imfx^t^A^)] of the response function Xa^a can be measured from the amount of energy 
(or momentum) transferred to the system for a given excitation sequence of the form 
(JHJ). For a quasi-monochromatic excitation at frequency Au, the transferred energy is 
indeed proportional to _E a b s oc Aulm\xAiA(Au)]. For a system at thermal equilibrium 
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at a temperature T, the fluctuation-dissipation theorem jT7j then relates Im[x j4 t^(Ao;)] 
to the corresponding dynamic structure factor S a ^ a (Auj) — fdt (A^(t) ^4(0)) e lAuJt : 

ImbuUM] = S A , A (Au;) (l - e~ hA ^ T ) . (12) 
In our case, S A t A can be written as: 

S AU {Au) = f n cos 2 9 S n (q, Alu) + f 2 s sin 2 9 S s (q, Au), (13) 
in terms of the density and spin dynamical structure factors S njS (q, Au) defined as: 

S „( q ,A.)=/ dtdr ( fl(r , t ) fl (0,0))e~0 (14) 

S s (q,Aw) = Jdtdr(s z (r,t)S z (Q,Q))^ l{cir ~ AuJt) - (15) 

The contributions of respectively the total and the spin density can be then 
isolated by measuring the response function Im[x^t^] for different values of the angle 9 
between the two polarizations. Knowledge of the total and spin density structure factors 
immediately gives the same- and opposite-spin density correlations: 



S\\ )± (q,Auj) = 




Differently from the dynamic structure factor, the low-temperature response function 
Xa^a is generally weakly dependent on T, and is well approximated by its T = value. 
Under this approximation, by setting T = in ()12j) . it is immediate to see that a 
measure of Xaia provides a good approximation of the most relevant zero temperature 
value of Sj^a(Auj). The static structure factor S A f A = (A^tyAit)) is obtained from the 
corresponding dynamical one S A ^ A (Au) by frequency integration: 

Sn. s (q) = S n , s (q, Aw) e-^K (17) 

J 27T 

4. Spontaneous scattering 

So far, we have considered the case of stimulated Bragg scattering, where two laser beams 
were sent on the cloud and scattering out of the first was stimulated by the presence 
of the second one. Now we shall analyse a different configuration, where information 
on the spin and density correlations in the atomic gas can be directly retrieved without 
invoking the fluctuation-dissipation theorem. Let us consider a spontaneous Bragg 
scattering process in which a single, linearly polarized, coherent laser beam is sent onto 
the atomic cloud with a wavevector ki and amplitude Ei = E° e\ and the intensity of 
spontaneously scattered light with wavevector k 2 , linear polarization e 2 , and frequency 
u>2 = ck 2 — Ui — Au is measured as a function of the scattering wavevector q = k 2 — kx 
and the polarizations ei )2 . 

For the sake of simplicity, we start from the case where no external field is present, 
so that the light-matter coupling Hamiltonian (JBJ) can be used. If the laser frequency uj\ 
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is far enough from the resonance frequency u ge , the atomic cloud is optically dilute and 
optically thin. Multiple scattering processes can then be neglected. Within this Born 
approximation, the (vector, operator-valued) scattering amplitude can be written as: 

E(k 2 ) ~ ^-V h2 (f n n(c L ,Aco)e 1 + if s e 1 x S(q,Aw)), (18) 

where k 2 is the wavevector of the scattered photon, n(q) and S(q) are the Fourier 
components of the total density and the spin density, and the projector Vk 2 projects 
onto the plane orthogonal to k 2 in order to account for the transverse nature of the e.m. 
field. This expression has a transparent physical intepretation in terms of a fluctuating 
dielectric susceptibility of the atomic gas jlj. The total density corresponds to the 
scalar, isotropic, component of the dielectric tensor, while the spin density gives an 
optical activity which rotates the polarization plane of the light. Fluctuations at a 
wavevector q of the susceptibility tensor scatter light from the incident mode at ki into 
the one at k 2 . Note that no nematic term of the form discussed in j3] for spin-1 atoms 
can be found in the present spin-1/2 case because of symmetry arguments. 

As we are considering atomic clouds whose density profile is smooth, the spatial 
Fourier transform of the average total and spin densities are non-vanishing only at 
wavevectors q, much smaller than the inverse cloud size 1/R. This implies that the 
coherent scattering amplitude (E(k 2 )) is not vanishing only in a small cone of aperture 
X/R around the incident beam direction, corresponding to the standard diffraction on 
a refractive body of size R. In the following, we shall place ourselves well outside this 
cone, where only incoherent scattering occurs. 

Interesting information on the density fluctuations in the cloud can be retrieved 
from the properties of the incoherently scattered light, in particular from the dependence 
of its intensity on the polarization vectors ei j2 . In practice, the scattered intensity J sc 
into the polarization state e 2 can be isolated by making the scattered light at k 2 to pass 
through a linear polarizer before reaching the detector. As e 2 is by definition orthogonal 
to k 2 , the intensity of the scattered light can be written as: 

J sc = ((e 2 -Et(k 2 ))(e 2 -E(k 2 )))= (19) 

= [f 2 n cos 2 9 S n (q, Au) + f 2 s sin 2 9 S s (q, Aw)] I mc (20) 

where the quantization axis z for the spin has been again chosen along the e 2 x ei 
direction, L mc = \E°\ 2 is the incident intensity and 9 is the angle between the linear 
polarization vectors ei )2 . Note the strict analogy between the result (j2*U|) for spontaneous 
scattering, and formula (fT3|) for the stimulated one. Although different quantities are 
experimentally measured in the two cases (the response function in the stimulated 
scattering, the structure factor in the spontaneous one), the atomic operator involved 
in the scattering process is the same, as one can also verify by comparing (J2J) to the 
expression for e 2 • E(k 2 ) obtained from (fTHj). As happens in the stimulated scattering 
case, the contribution of the density and spin structure factors to (J20|) can be isolated 
simply by repeating the scattering experiment for different values of the angle 9 between 
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the polarization vector of the incident beam and the polarization direction of the 
polarizer used for analysing the scattered light. 

In order to directly access the static structure factor S n>s (q), one has to detect 
the whole spectrum of the scattered light (as in X-ray scattering from crystals [TT?] ) 
or, better, to perform the experiment using a short pulse containing a wide range of 
frequencies as the incident beam. With respect to the stimulated Bragg scattering 
discussed in the previous section, the present technique has therefore the advantage 
of not requiring the repeated measurement for different detunings Au; in order to get 
the static structure factor f. This is easily understood in physical terms: a scattering 
experiment performed with a short pulse compared to the characteristic time-scale for 
the atomic motion provides an instantaneous measurement of the atomic positions, i.e. 
of the static structure factor. 

A major disadvantage of the spontaneous Bragg scattering technique is the 
limitation on the maximum number of photons (on the order of the number of atoms in 
the cloud) that can actually be scattered off the cloud before information on the internal 
state of the atoms is destroyed by the spontaneous emission processes. 

5. Role of external magnetic fields 

As most experiments on superfluidity in two-spin component Fermi atoms require 
an external magnetic field Bq to enhance atom-atom interactions via the Feshbach 
resonance effect [21 El, it is important to extend the discussion to this experimentally 
relevant case. The geometrical issues are not as simple as discussed in secEJ so we shall 
limit ourselves to a few specific examples of stimulated Bragg scattering processes and 
we shall postpone a more complete analysis to a forthcoming publication. 

Consider first where the magnetic field is weak enough not to mix the 

F g = 1/2 ground state with other hyperfine components; this approximation is expected 
to be valid provided the Zeeman energy [iBq (/i is the magnetic moment of the atom 
in the ground state) is much smaller than the hyperfine splitting A^f- For the sake of 
simplicity, we assume that both Bragg beams are polarized on the plane orthogonal to 
the magnetic field direction (which is taken as the z axis). In the circular ± basis, the 
Hamiltonian Q is replaced by: 



| Note that the stimulated Bragg scattering experiment discussed in secdl can not be used in the 
short pulse regime in order to directly measure the static structure factor. An extra Aw factor is in 
fact present in the proportionality law between the absorbed energy £" a bs and the quantity of interest 
Im[x J 4tyi(Aa;)], so that the total absorbed energy in a pulsed experiment is trivially determined by the 
/-sum rule [T%] : 





ij=± 



(22) 




(21) 



rather than by the static structure factor. 
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= \j dr { [(/++ + /+-)*+(*) + (/-+ + f-)L(r)]h(r) + (23) 
+ [((/++ " U-)h(r) + (/-+ - /-)/-(r)] &(r)} (24) 

where i±(r) = E±(r) E±(r) are the local light intensity operators in the a± polarization 
state. The oscillator strenghts for the j = ± Zeeman state interacting with 
o~i=± polarized light have a non-trivial dependence on the light frequency ul and the 
resonance frequencies of the Zeeman-split transitions. For the case of a transition 
to a F e = 1/2 state, no o~± transition can be driven on atoms in respectively the ± 

state, so / ++ = / = 0. On the other hand, /+_ = \d\ 2 /h(uL — w e ,+i/2 + Wg,-) an d 

/_+ = \d\ 2 /h{ujL — w e ,-i/2 + u g,+)- F° r the more complicate case of a transition to a 
F e = 3/2 state, one has to take into account the relevant Clebsch-Gordan coefficients, 
so that /++ = \d\ 2 /h(u L - ^ e , +3 /2 + ^g,+), /-+ = \d\ 2 /3h(u L - Ld e _ 1/2 + u 9t+ ) and 

analogous formulas hold for the / H and / oscillator strengths starting from the g, — 

initial state. Here, 0J g ,i=± and w e ,i=+3/2,+i/2,-i/2,-3/2 are the energies of the Zeeman-split 
sublevels of the ground and excited states. 

The analysis of sec|Hl can then be easily repeated for this Hamiltonian: the main 
difference is in the perturbation operator which now reads: 



+ + (/-+ + /->" A(q)+ (25) 



i = I 
2 

+ \ [(/++ " U-)e' 16 + (/-+ - f-)e ie ] &(q). (26) 

As previously, 9 is the angle between the linear polarizations of the two beams. A 
non- vanishing coupling to the spin is present as soon as either / ++ ^ / + _ or /_ + ^ f 

Of course, this framework is directly generalized to the case where the two Zeeman 
sublevels of the ground state are replaced by a pair of sublevels belonging to a higher 
angular momentum state, as in recent experiments with K atoms trapped in their 
F = 9/2, mp = —9/2,-7/2 states j^J- The only difference is in the Clebsch-Gordan 
coefficients appearing in the expression of the oscillator strengths fy. In 40 K atoms, the 
hyperfine splitting Ahf ~ 870 MHz is in fact much larger than the Zeeman splitting at 
the B Q ~ 200 G Feshbach resonance nB w 280 MHz Ahf- Further extension to the 
case where the two atomic states belong to different hyperfine states is straightforward; 
this was indeed the case of the state-selective imaging scheme used to study vortices in 
two-component Bose-Einstein condensates 

The case of 6 Li is more complicated: around the B ~ 850 G Feshbach resonance 
currently used in experiments on fermionic superfluidity, the Zeeman splitting u,Bq ~ 
1.2 GHz is substantially larger than the ground state hyperfine splitting ~ 225 MHz 
which makes the level structure significantly different from the one discussed upto now. 
We postpone a complete discussion of these issues to further investigation, and we 
here limit ourselves to some general considerations. In this regime, the atomic internal 
eigenstates are well approximated by the ones of the Zeeman Hamiltonian, and the 
hyperfine coupling is only a small perturbation. In particular, atoms are trapped in 
the two lowest-lying ± states in which the electron is always in the total spin state 
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J = 1/2, mj = —1/2 and the nucleus is in respectively the 1 = 1, mj = 1,0 spin 
state • The fact that the electron is almost decoupled from the nucleus and is in the 

J = 1/2, mj = —1/2 state implies that the oscillator strenghts / and /_+ for o~_- 

polarized light are very weak (the matrix elements are a factor of the order n = A jjp / [iB Q 
smaller). The coupling to the spin is mainly due to the difference Af = f ++ — /+_, 
and is provided by the hyperfine coupling in both the ground and the excited states. 
This in fact splits the resonance frequencies for the o~ + transitions starting respectively 
from the ± states. To make Af the largest, one has to choose a laser frequency close 
to resonance, but a compromise has to be made with spontaneous emission. 

6. Conclusions 

In summary, we have shown how it is possible to take advantage of the light polarization 
degrees of freedom to separately address the total density and the spin density of a 
two-spin component atomic gas. This can be exploited in polarization-sensitive Bragg 
scattering techniques so as to measure both the density and the spin structure factors 
of the gas and then to characterize the microscopic structure of the gas in its different 
phases. Stimulated Bragg scattering probes the response function of the gas to the 
optical potential created by the interference pattern of the two laser beams: one can act 
on either the total or the spin density simply by changing the light polarizations. The 
fluctuation-dissipation theorem then relates the absorbed energy to the structure factor 
of our interest. On the other hand, the density and spin structure factors can be directly 
probed in a spontaneous Bragg scattering experiment: the intensity of scattered light 
into a given polarization state is in fact proportional to a combination of the density 
and the spin structure factors, whose weights depend on the relative angle between the 
incident and scattered beam polarizations. 
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continuous stimulating exchanges. Stimulating discussions with R. Hulet and C. Regal 
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